
> # Read in the data
> load("REStatSubmission_FRB.RData")

> # Only retain the observations that will be displayed (Census Tracts with relative incomes
> # of between 60 and 140 percent)
> myData <- subset(myD .... [TRUNCATED] 

> # Create a smaller data frame that contains only the variables that will be plotted
> graphData <- data.frame( count=myData$COUNT,
+                 .... [TRUNCATED] 

> #######################################################################
> ###  Figure 1 (This is the only figure in the paper)
> ################### .... [TRUNCATED] 

> # Set up display as a matrix of 6 graphs (3 x 2) and set margins around the graphs
> par(mfrow=c(3,2),
+     mar=c(4,3,2,1) )

> # This function calculates weighted means by specified group.  The groups used in the 
> # graphs will be the levels of relative income between 60 a .... [TRUNCATED] 

> # This function demeans a variable (x), using a weighted mean (with weight equal to w),
> # by group.  The group used in the graphs will be the rela .... [TRUNCATED] 

> # Demean the six variables that will be plotted
> graphData$delrate.dm <- demean.variable( x=graphData$delrate,
+                                    .... [TRUNCATED] 

> graphData$hidti_20046.dm <- demean.variable( x=graphData$hidti_20046,
+                                              w=graphData$loans_20046,
+      .... [TRUNCATED] 

> graphData$hirate_20046.dm <- demean.variable( x=graphData$hirate_20046,
+                                               w=graphData$loans_20046,
+   .... [TRUNCATED] 

> graphData$gse_20046.dm <- demean.variable( x=graphData$gse_20046,
+                                            w=graphData$loans_20046,
+            .... [TRUNCATED] 

> graphData$shr.dep.in_20046.dm <- demean.variable( x=graphData$shr.dep.in_20046,
+                                                   w=graphData$loan .... [TRUNCATED] 

> graphData$pur.dep.in_20046.dm <- demean.variable( x=graphData$pur.dep.in_20046,
+                                                   w=graphData$loan .... [TRUNCATED] 

> #############################################################
> ### PANEL A - Delinquency Rate
> ################################################### .... [TRUNCATED] 

> temp <- subset( data.frame(temp), 
+                 group.id>60 & group.id<=140 )

> plot(x=temp[,1],
+      y=temp[,2],
+      xlim=c(60,140),
+      #ylim=c( min(temp[,2])*0.9, max(temp[,2]*1.1) ),
+      type="o",
+      pch=20,
+ .... [TRUNCATED] 

> legend("topright",
+        c("CRA Threshold","GSE Threshold"),
+        lty=c(3,2),
+        cex=1,
+        xjust=1,
+        yjust=0.5)

> axis(2,las=1)

> axis(1,las=1,at=c(60,70,80,90,100,110,120,130,140))

> abline(v=90,lty=2)

> abline(v=80,lty=3)

> #############################################################
> ### PANEL B - High Payment-to-Income Share
> ####################################### .... [TRUNCATED] 

> temp <- subset( data.frame(temp), 
+                 group.id>60 & group.id<=140 )

> plot(x=temp[,1],
+      y=temp[,2],
+      xlim=c(60,140),
+      type="o",
+      pch=20,
+      main="(b) % High PTI",
+      bty="l",
+      ylab .... [TRUNCATED] 

> axis(2,las=1)

> axis(1,las=1,at=c(60,70,80,90,100,110,120,130,140))

> abline(v=90,lty=2)

> abline(v=80,lty=3)

> #############################################################
> ### PANEL C - Share of Higher-Priced Loans
> ####################################### .... [TRUNCATED] 

> temp <- subset( data.frame(temp), 
+                 group.id>60 & group.id<=140 )

> plot(x=temp[,1],
+      y=temp[,2],
+      xlim=c(60,140),
+      type="o",
+      pch=20,
+      main="(c) % Higher Priced",
+      bty="l",
+      .... [TRUNCATED] 

> axis(2,las=1)

> axis(1,las=1,at=c(60,70,80,90,100,110,120,130,140))

> abline(v=90,lty=2)

> abline(v=80,lty=3)

> #############################################################
> ### PANEL D - GSE Sales (2004-2006)
> ############################################## .... [TRUNCATED] 

> temp <- subset( data.frame(temp), 
+                 group.id>60 & group.id<=140 )

> plot(x=temp[,1],
+      y=temp[,2],
+      xlim=c(60,140),
+      type="o",
+      pch=20,
+      main="(d) GSE Sales",
+      bty="l",
+      ylab= .... [TRUNCATED] 

> axis(2,las=1)

> axis(1,las=1,at=c(60,70,80,90,100,110,120,130,140))

> abline(v=90,lty=2)

> abline(v=80,lty=3)

> ###########################################################################
> ### PANEL E - Lending by Depositiories Within Assessment Areas (2004-2 .... [TRUNCATED] 

> temp <- subset( data.frame(temp), 
+                 group.id>60 & group.id<=140 )

> plot(x=temp[,1],
+      y=temp[,2],
+      xlim=c(60,140),
+      type="o",
+      pch=20,
+      main="(e) Lending: Depository - In",
+      bty="l ..." ... [TRUNCATED] 

> axis(2,las=1)

> axis(1,las=1,at=c(60,70,80,90,100,110,120,130,140))

> abline(v=90,lty=2)

> abline(v=80,lty=3)

> #####################################################################################
> ### PANEL F - Share of Purchases by Depositories within Asse .... [TRUNCATED] 

> temp <- subset( data.frame(temp), 
+                 group.id>60 & group.id<=140 )

> plot(x=temp[,1],
+      y=temp[,2],
+      xlim=c(60,140),
+      type="o",
+      pch=20,
+      main="(f) Purchases: Depository - In",
+      bty= .... [TRUNCATED] 

> axis(2,las=1)

> axis(1,las=1,at=c(60,70,80,90,100,110,120,130,140))

> abline(v=90,lty=2)

> abline(v=80,lty=3)

> dev.off()
null device 
          1 

> sink()  # Close log file
